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LONG-TERM  GOAL 

The  long-term  goal  of  this  work  is  to  develop  an  accurate  and  efficient  LES  (large-eddy  simulation) 
model  of  turbulent  flows  in  oceanic  boundary  layers. 

OBJECTIVES 

The  major  objectives  to  be  achieved  are: 

(1)  Acquire  better  understanding  of  specific  properties  of  oceanic  boundary  layers,  particular  attention 
being  paid  to  statistics  of  small-scale  coherent  structures  and  their  adequate  representation  by  subgrid- 
scale  (SGS)  models. 

(2)  Implement  the  Smagorinsky  and  dynamic  SGS  models  and  validate  them  for  oceanic  boundary 
layers  through  comparison  with  direct  numerical  simulations  (DNS)  and  autonomous  underwater 
vehicle  (AUV)  field  data. 

(3)  Develop  and  test  a  SGS  model  capable  of  accurate  parameterization  of  the  effects  of  small-scale 
coherent  structures. 

APPROACH 

The  work  involves  theoretical  analysis,  numerical  computations,  and  comparison  with  field 
measurements.  The  primary  research  tools  are  three-dimensional  DNS  and  LES  models,  which  allow 
simulation  of  turbulent  layers  in  shallow  water. 

WORK  COMPLETED 

The  three-dimensional  numerical  model  developed  by  Slinn  &  Riley  (1998)  has  been  modified  to 
describe  turbulent  penetrative  convection  in  shallow  water  caused  by  the  passage  of  a  cold  air  front. 
This  flow  was  chosen  because  of  an  extensive  database  of  field  measurements  collected  during  AUV 
experiments  in  December  1998.  We  started  with  two  high  resolution  Direct  Numerical  Simulations 
performed  at  Rayleigh  number  Ra=7  xlO*  with  different  background  density  stratification  (Richardson 
number  Ri  =  2  and  10).  Further,  we  have  implemented  the  standard  LES  Smagorinsky  model  and  the 
scaling  version  of  the  dynamic  LES  SGS  model  proposed  for  turbulent  convection  by  Wong  &  Lilly 
(1994).  The  models  have  been  tested  and  validated  through  comparison  with  results  of  DNS.  In  the 
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course  of  the  work,  we  have  found  that  the  high  Rayleigh  numbers  observed  in  typical  oceanic  flows 
can  be  achieved  only  upon  further  modification  of  the  model,  in  partieular,  use  of  approximate 
boundary  eonditions  at  the  upper  surfaee  beeomes  expedient. 

RESULTS 

We  eonsider  the  flow  in  a  horizontal  layer  eonfined  between  a  smooth  no-slip  bottom  and  rigid  firee- 
slip  surfaee  at  the  top  (see  Figure  1).  The  flow  is  assumed  homogeneous  in  both  horizontal  direetions 
so  that  periodie  boundary  eonditions  ean  be  applied.  The  flow  is  generated  by  the  mechanism  of 
thermal  eonvection,  whieh  arises  due  to  heat  exehange  (eooling)  at  the  upper  surfaee.  Stable 
baekground  stratifieation  is  present  initially,  eaused  by  a  baekground  temperature  gradient. 


1.  Schematics  of  the  model  problem. 


The  flow  evolution  is  illustrated  in  Figure  2.  Instantaneous  density  and  veloeity  fields  ealculated  in  the 
DNS  run  with  the  stronger  baekground  stratifieation  are  shown  in  a  vertieal  eross-seetion  from  the 
three-dimensional  simulation.  It  ean  be  seen  how  the  initial  diffusive  layer  (Figure  2a)  beeomes 
unstable  (Figure  2b)  and  ereates  turbulent  eonveetive  motion  (Figure  2c).  The  turbulenee  penetrates 
the  layer  and,  after  some  time,  reaehes  the  bottom  boundary,  forming  a  turbulent  layer  aeross  the  water 
eolumn  (Figure  2d). 

Our  caleulations  have  eonfirmed  the  eonclusion  made  previously  for  Rayleigh-Benard  eonveetion  (see, 
e.g.,  Edison,  1985)  that  the  aeeuraey  of  simulations  depends  erueially  on  the  numerieal  resolution  of 
the  diffusive  boundary  layer  at  the  upper  surfaee.  In  the  DNS  runs,  the  vertieal  grid  size  was  deereased 
almost  tenfold  in  this  region.  In  LES  eomputations,  the  vertical  grid  size  ean  be  inereased  in 
eomparison  to  DNS  only  upon  replacing  the  boundary  eondition  of  a  given  heat  flux  at  the  surfaee  by 
its  approximation  over  the  grid  eontrol  volume  elosest  to  the  boundary.  In  the  first  EES  simulations 
presented  in  this  report,  we  retain  the  exact  boundary  condition  and  apply  the  SGS  model  in  the 
horizontal  directions. 


t=0.83 


2.  Evolution  of  the  density  field  in  a  two-dimensional  cross-section  from  the 
three-dimensional  DNS  experiments  at  different  times. 


Two  SGS  models  have  been  implemented  in  the  eode.  One,  the  so-ealled  scaling  formulation  of  the 
dynamic  (Germano)  model  proposed  for  turbulent  convection  by  Wong  &  Lilly  (1994),  has  been  found 
to  produce  relatively  weak  SGS  stresses.  By  comparison,  the  classical  Smagorinsky  model  with  fixed 
turbulent  Prandtl  number,  Prx=l/2,  has  proven  robust  and  capable  of  predicting  the  average  flow 
characteristics  in  close  agreement  with  DNS.  Comparison  between  LES  and  DNS  is  illustrated  in 
Figure  3.  Figure  4  shows  the  time  evolution  of  the  horizontally  averaged  turbulent  kinetic  energy  as  a 
function  of  height  from  the  boundaries. 


3.  Comparison  of  TKE  and  potential  energy  contributions  between 
LES  and  DNS  simulations. 


IMPACT/APPLICATIONS 

An  important  potential  benefit  of  our  work  for  soeiety  is  that  it  serves  to  develop  and  validate  an 
aeeurate  and  effieient  LES  model  for  turbulent  flows  in  shallow  water.  The  model,  ineluding  effeets  of 
small-scale  coherent  structures  on  SGS  stresses,  will  enhance  our  ability  to  simulate  such  important 
processes  as  diapycnal  mixing  in  shallow  water,  wave  boundary  layers,  or  beach  processes. 
Furthermore,  our  work  contributes  to  better  understanding  of  properties  of  turbulent  flows  in  the 
shallow  ocean,  in  particular,  turbulent  penetrative  convection  caused  by  surface  cooling. 


4.  TKE  distribution,  averaged  in  x-y  plane  from  a  DNS  computation, 
TRANSITIONS 


One  important  transition  in  our  work  has  been  eaused  by  the  faet  that  we  extensively  used  the  field 
observations  from  the  adverse  weather  experiments  eonducted  in  a  related  projeet  to  validate  the 
numerical  models.  Accordingly,  the  simulations  have  focused  on  turbulent  convection  generated  by 
surface  cooling  rather  than  bottom  boundary  layers  with  prescribed  free-stream  conditions  as  was 
initially  planned  for  the  project. 

RELATED  PROJECTS 


The  work  is  carried  out  in  conjunction  with  field  measurements  under  grant  N00014-96-1-5023. 
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